ADP-12-10/T777 



Structure and Coulomb dissociation of 23 O within the quark-meson 

coupling model 

R. Chatterjee 1 , R. Shyam 2 ' 3 , K. Tsushima 24 , and A. W. Thomas 2 
1 Physics Department, Indian Institute of Technology, Roorkee, India 
1 Centre for the Subatomic Structure of Matter (CSSM), 
School of Chemistry and Physics, University of Adelaide, SA 5005, Australia 
3 Saha Institute of Nuclear Physics, 1/AF Bidhan Nagar, Kolkata, India and 
4 International Institute of Physics, Federal University of Rio Grande do Norte, Natal, Brazil 

(Dated: March 2, 2013) 

Abstract 

We study the ground-state structure of the near drip-line nucleus 23 within the latest version of the 
quark-meson coupling model. It predicts a 22 O(0 + )®n2si/2 configuration for the ground state of 23 0, with 
a one neutron separation energy in close agreement with the corresponding experimental value. The wave 
function describing the valence neutron-core relative motion obtained within this model, has been used to 
calculate the Coulomb dissociation of 23 O on a lead target at a beam energy of 422 MeV/nucleon. The 
experimental neutron-core relative energy spectrum and the total one-neutron removal cross sections are 
well described by our calculations. The widths of the longitudinal momentum distributions of the 22 
fragment are found to be broad, which do not support the formation of a neutron halo in this nucleus. 

PACS numbers: 25.60.Gc, 12.39.Ki, 24.85.+p 



1 



In recent years a lot of effort has been devoted to study the structure of the neutron rich 23 
nucleus that has 15 neutrons (AT). This is driven by several unusual observations. One of them 
is that the neutron drip-line for the oxygen isotope is located at 24 (see e.g. Refs. [1]), and not 
at 28 0, the doubly magic nucleus that has been found to be unbound. This suggests that a new 
magic number appears at N = 16 in neutron rich nuclei. The reason for this observation has been 
investigated in several shell model based studies [2 j. For example, Otsuka et al. M\ have suggested 
that it should be attributed to the spin-isospin dependent part of the nucleon-nucleon force yielding 
a strongly attractive interaction between a proton in 1 J 5 / 2 orbit and a neutron in 1 J3/2 orbit, which 
in turn increases the gap between 2s 1 / 2 and U3/2 orbitals for the N=16 oxygen isotope. 

As the neutron drip-line is approached the nuclei experience weakening of the neutron binding 
energies, which leads to some special effects. The sudden rise of the interaction cross sections 
found in Refs. |3|] for several lighter isotopes closer to the neutron drip-line was attributed to the 
extended density distribution(s) of the valence neutron(s) usually referred to as a neutron halo Q4J] . 
For oxygen isotopes a sudden rise of about 12% was observed D5J in the interaction cross sections 
between the 22 and 23 0. This led these authors to suggest 23 as a candidate for a neutron-halo 
system. However, there are several other observations that are not consistent with this picture. 
The valence neutron in 23 is relatively strongly bound, with a one neutron separation energy 
(S n ) of 2.74 MeV, which is in contrast to the well established one neutron halo nuclei like, n Be 
and 19 C, that have S„ of only about 0.5 MeV. Thus the last neutron in 23 should not extend to 
larger distances. Furthermore, the longitudinal momentum (py) (with respect to the reaction plane) 
distributions (LMD) of the 22 O fragment have been measured [6-I2] in the breakup reactions of 
23 on lighter target nuclei, where the widths of these distributions are found to be in excess of 100 
MeV/c. This is about 2-3 times larger than those observed in the cases of 11 Be and 19 C induced 
reactions (see, e.g. Ref. [10]). A narrow LMD is a direc t sig nature of the development of the halo 
structure (see, e.g. the discussions presented in reviews [UjJ]). 

Moreover, there has been a disagreement in the assignment of the ground state spin-parity (J*) 
of 23 nucleus among various authors. In Ref. Q6Q, from the LMD measurements of the 22 and 
21 fragments in the 23 induced reaction on a carbon target, it was concluded that the data support 
a J" of | + . This conclusion, however, was not supported by the measurements of the 23 breakup 
reactions reported in Ref. Q8 , 



1211 where a J n of ^ + was found to be consistent with the data. In 



Ref. [13] the interaction cross sections of 22 and 23 isotopes were measured on a carbon target 
and it was concluded that these data are more in agreement with a 22 O + valence neutron in the 



2*1/2 picture of the 23 nucleus. This would imply a J" of j + . These conflicting experimental 
results call for a more rigorous theoretical investigation of the 23 ground state configuration. 
Furthermore, there is a need to analyze the data on the breakup reactions of 23 within a more 
microscopic theoretical model than what has been done so far. 



The aim of this Rapid Communication is to investigate the ground state structure of the 23 
nucleus within the quark-meson coupling (QMC) model, which is a quark-based model for nuclear 
matter, finite nuclei and hypernuclei [|14l-ll9|l . In the QMC model, quarks in the non-overlapping 
MIT bags interact self consistently with isoscalar-scalar (cr) and isoscalar-vector (to) mesons in the 
mean field approximation. The explicit treatment of the nucleon internal structure is a key point 
of this model and it represents an important departure from quantum hadrodynamics (QHD) [20]. 
The self-consistent response of the bound light quarks to the mean cr field leads to a new saturation 
mechanism for the nuclear matter II 1411 . Although formulated as a relativistic mean field theory 
with just a few parameters (determined from the properties of nuclear matter), QMC has been 
shown to lead to a remarkably realistic Skyrme force [21]. 



In the present calculation, we use the latest version of the QMC model (to be referred to as 
QMC-III) |Q] that includes, in particular, the self-consistent effect of the mean scalar field on the 
familiar one-gluon exchange hyperfine interaction [|22ll that in free space leads to the N - A and 
Z-A mass splitting. With this development, QMC model has been able to explain the properties of 
A hypernuclei for the s-state rather well, although the p- and J-states tend to be underbound. It also 
maintains the very natural explanation of the small spin-orbit force in the A-nucleus interaction 



that was found in earlier version of the QMC model (to be referred to as QMC-I) ||16L 1 1 VL 



230. In 



QMC-III, while the quality of results for A and S hypernuclei is comparable to that obtained in 



-I Hi 



QMC-I 111711 . no bound states for the Z states 111911 were found in medium and heavy mass nuclei. 
This finding, which is a consequence of the extra repulsion associated with the increased one- 
gluon-exchange hyperfine interaction in medium, is in agreement with the non-observation of such 
states experimentally. The QMC-III model has recently been used [24j] to study the production of 
S~ hypernuclei via the (K~,K + ) reaction, which is currently of great interest at the J-PARC facility 
in Japan. 



In order to calculate the properties of the finite nuclei, we construct a simple, relativistic shell 
model with self-consistent scalar and vector mean fields. The Lagrangian density for a nuclear 
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system in the QMC model is written as 111611 : 



Lqmc = ^ N {r)\iy ■ d - M N {cr) - (g a (o(r) + g p ^-b(r) + U\ + r^)A(r) )y ]tff N (r) 
-^[(Vcr(r)) 2 + mlcr(r) 2 ] + ^[(Vco(r)) 2 + m>(r) 2 ] 
+^[(VKr)) 2 + m 2 p b(r) 2 ] + ^(VA(r)) 2 . 



(1) 



Here ^(r), £(r), <y(r) and A(r) are, respectively, the nucleon, the p meson, the co meson and 
Coulomb fields, while m a , m w and m p are the masses of the cr, co and p mesons. g w and g p are 
the co-N and p-N coupling constants that are related to the corresponding (u,d)-quark-a», gl, and 
(w, rf)-quark-p, g q p , coupling constants as g^ = 3gtj and g p = g p , and e is the proton charge. 
The following equations of motion are obtained from the Lagrangian density Eqs. ©: 

[iy ■ d - M N (cr) - ( g^{r) + g p ^j-b(r) + (1 + if )A(r) ) 7o ]^(r) = 0, (2) 



(-V 2 + m 2 )o-(r) = gcrC N (p-)p s (r), 



(3) 



(-V 2 +m>(r) = ££p 3 (r), 



(4) 



(-V 2 )A(r) = e Pp (r), 



(5) 



where, p,.(r), p B (r) = p p (r) + p„(r), p 3 (r) = p p (r) - p„(r) andp p (r) (p„(r)) are the scalar, baryon, 
third component of isovector, and proton (neutron) densities II 1711 . On the right hand side of Eq. ©, 
a new and characteristic feature of QMC appears, arising from the internal structure of the nucleon, 
namely, go-C N (o-) = - dM ^ cr) , where go- = go-(o~ = 0). We use the density dependent nucleon mass 
M N (cr) as parameterized in Ref. IU9I1 . 

The coupled non-linear differential Eqs. (2)-(5) can be solved by an standard iteration procedure 
as discussed in, e.g., Refs. j25\ \. The coupling constants go-, g u , g p , and masses m^, m w and m p 
have been taken to be the same as those given in Ref. [19Q. To calculate the nuclear levels we use a 



relativistic she 



we refer to [116 



1 model. For a more detailed description of the finite nuclei within the QMC model 



ISO- 



In our calculations for the 23 nucleus, the configuration where valence neutron lies in the 2si/2 
level (with 22 remaining as inert core), yields a one neutron binding energy of 2.86 MeV. This 
is very close to the experimental value of 2.74 MeV for this quantity Q26[|. We have also tried a 



scheme to get a valence neutron spin-parity of | by filling the last neutron orbit (2si/ 2 ) with two 
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FIG. 1. Neutron (solid line) and proton (dashed line) density distributions for 23 O in QMC-III model. 



neutrons. This, however, makes the 23 isotope unstable in our model. Therefore, we confirm that 
the ground state of this nucleus is consistent with a 22 O(0 + )®2.yi/2 configuration with a one neutron 
separation energy of 2.86 MeV. It should be mentioned here that using the QMC-I model, we find 
a one neutron binding energy of 4.34 MeV within the same ground state configuration. Thus, the 
extra effects put into the QMC-III model (e.g., one-gluon exchange hyperfine interaction) do seem 
to be important to describe self-consistently the structure of the drip-line nuclei. This point will 
be investigated in more detail in a future study. 

In Fig. 1, we show the neutron and proton density distributions for 23 obtained within the 
QMC-III model. We see that the central density of the neutron distribution is clearly larger than 
that of the proton. At the same time the surface diffuseness of the two distributions are almost 



identical. This is consistent with the criteria for a neutron skin [|27L 



280 . The neutron skin thickness, 

which is the difference between the neutron and proton root mean square radii (RMS) 12911 is found 



to be 0.46 fm, which is in agreement with the value reported for this nucleus in Ref. [12811 where 
calculations have been performed within a QHD type of relativistic mean field model. In neutron 
rich isotopes of several medium to heavy mass nuclei, a neutron skin of similar thickness has been 
reported by several authors (see, e.g. Refs. BOIl ). 

In Figs. 2a and 2b, we show the valence neutron-core potential and the corresponding wave 
function, respectively, for the 23 nucleus calculated within the QMC-III model. The potential is 
the sum of scalar and vector fields and the wave function is only the upper component of the spinor, 
as the corresponding lower component is at least two orders of magnitude smaller. For comparison, 
we also show a Woods-Saxon (WS) potential and the respective wave function for this system. 
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FIG. 2. The QMC-III (solid line) and Woods-Saxon (dashed line) valence neutron-core potentials [part (a)] 
and corresponding wave functions [part (b)] for the 23 nucleus. 



This potential has been parameterized by adjusting its depth to reproduce the experimental value 
of S„ with the radius and diffuseness parameters of 1.15 fm and 0.5 fm, respectively. We note that 
while the QMC-III potential is somewhat deeper than the WS potential at smaller distances, they 
are similar in the surface region. However, the wave functions generated by the two potentials are 
almost identical at all radii. 

The wave functions calculated above have been used to investigate the Coulomb dissociation 
(CD) of the 23 nucleus on Pb target, for which some data exists [12] on the valence neutron-core 
relative energy differential cross section and the total one-neutron removal cross sections. The CD 
process has a number of advantages. It is free from the uncertainties associated with the nuclear 
interactions. Moreover the inelastic breakup mode (also known as stripping or breakup fusion) 
is absent in the pure Coulomb dissociation reactions. Therefore, this process is ideally suited for 
probing the structure of the projectile. 

The model used to calculate the CD cross sections of the 23 O nucleus is described in detailed in 



Refs. p 1 M33H - Therefore, we give here only a brief sketch of it. This theory is formulated within 
the post-form finite range distorted wave Born approximation (FRDWBA), where the electromag- 
netic interaction between the fragments and the target nucleus is included to all orders and the 
breakup contributions from the entire nonresonant continuum corresponding to all the multipoles 



and the relative orbital angular momenta between the fragments are taken into account 113211 . Only 
the full ground state wave function of the projectile, of any orbital angular momentum configu- 
ration, enters in this theory as input. Thus, unlike most of the theoretical models described in a 



recent review of breakup theories V34n . this model does not require the knowledge of the positions 



and widths of the continuum states. Hence, it is free from the uncertainties associated with the 



multipole strength distributions [35] that may exist in some of these theories. 

We consider a breakup reaction (a + A —> b + c + A), where the projectile a breaks up into 
fragments b (charged) and c (uncharged) in the Coulomb field of a target A. The differential cross 
section for the relative energy distribution for this reaction is given by 

dE bc J^ A \^( 2l+ 1) I h6 



where v aA is the a— A relative velocity in the entrance channel, Q. bc an d Cl aA are solid angles, /j. bc 
and n aA are reduced masses, and p hc and p aA are appropriate linear momenta corresponding to the 
b-c and a— A systems, respectively. 

The reduced amplitude, /?/„,, in the post-form FRDWBA is given by 

fo m = (cxp(yk c - aK)\V bc \^)^ b -\k b )x i ;\6k c )\x^\k a )) , (7) 

where k b , k c are Jacobi wave vectors of fragments b and c, respectively, in the final channel of 
the reaction, k a is the wave vector of projectile a in the initial channel and V bc is the interaction 
between b and c. &™ is the ground state wave function of the projectile with relative orbital angular 
momentum state / and projection m. In the above, K is an effective local momentum associated 
with the core-target relative system, whose direction has been taken to be the same as the direction 
of the asymptotic momentum k b 131(1. a, 6, and y in Eq. (5), are mass factors relevant to the Jacobi 
coordinates of the three body system (see Fig. 1 of Ref. BUJ]). The wave functions and^ , 
are the distorted waves for the relative motion of b and c with respect to A and the cm. of the 
b-A system, respectively, with ingoing wave boundary condition and;ti +) (£ a ) is the distorted wave 
for the scattering of the cm. of the projectile a with respect to the target with outgoing boundary 
conditions. 

The first term in Eq. © contains the structure information about the projectile through the 
ground state wave function 0^ m , while the second term is associated with the dynamics of the re- 
action. For the pure Coulomb case, x b ~\k b ) and x { a\k a ) are replaced by the appropriate Coulomb 
distorted waves, and;f£ (dk c ) by a plane wave as the fragment c is uncharged (e.g. for a neutron). 
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FIG. 3. The calculated excitation energy spectra in the Coulomb dissociation of 23 O on a Pb target at 
422 MeV/nucleon beam energy obtained with QMC-III (solid line) and Woods-Saxon (dashed line) wave 



function. The data is taken from Ref. 11211 . The calculated cross sections are normalized to the peak of the 
experimental data in each case. 

This allows the second term of Eq. © to be evaluated analytically in terms of the bremsstrahlung 
integral 0611 . A more detailed description of the evaluation of the reduced amplitude (5 hn can be 



found in Refs. [31, 



37] 



It should be noted that in Eq. ©, the interaction between fragments b and c is treated to first 
order only. Because for relative energies of our interest there are no resonances in the n+ 22 
continuum, we expect this approximation to be valid. It is clearly a good approximation for the 
deuteron and the neutron halo systems [35]. For those cases where higher order effects of the 
fragment-fragment interaction are known to be nonnegligible, this model will have a limited ap- 
plicability. Nevertheless, we note that in calculating the relative energy spectra explicit information 
about the continuum strength distribution of the projectile is not required - the entire continuum is 
automatically included in this theory. 

In Fig. 3, the calculated energy differential cross section is shown as a function of the 
excitation energy (E x = Eb c + S n , where Et, c is the n- 22 relative energy). The results obtained 
with QMC-III (solid line) as well as a WS ground state (dashed line) wave functions are displayed. 
In getting these cross sections, an integration has been carried out over 6„ a up to the grazing angle 
of the reaction. In order to compare the calculations with the data of Ref. [12], the calculated cross 
sections need to be convoluted with the detector response function of the experiment. Since this 
is not available to us we have normalized our cross sections to the peak of the experimental jj- in 



^ 160 

I 120 



.Q 

E 



80 



a 

5 40 

ID 
T3 



QMC-III 




- - WS 






■ 






■ A 









19.4 19.5 19.6 19.7 

P„ (GeV/c) 



FIG. 4. The calculated longitudinal momentum distribution of 22 fragment in the Coulomb dissociation of 
23 on a Pb target at 422 MeV/nucleon beam energy obtained with QMC-III (solid line) and Woods-Saxon 
(dotted line) wave function. 

each case. This may be viewed as an alternative to the convolution procedure. 

We note that our calculations are able to describe the shapes of the relative energy spectrum 
quite well over the entire region of excitation energies. The sharp rise of the experimental cross 
section just after the threshold is quite well reproduced. Even the widths of the experimental 
distribution is well reproduced in our model. In contrast to this, the width of the distribution is 
grossly over predicted within a semiclassical breakup model using plane waves for the relative 
motion wave functions in the outgoing channel, as shown in Ref. [12]. Once optical potentials are 
used to describe the n + 22 relative motion, the situation improves. Nevertheless, given the well 
known uncertainty in this potential the agreement between theory and the data shown by these 
authors should be viewed with caution. 

The total electromagnetic one-neutron removal cross sections obtained in our CD model is 
79.26 mb and 77.15 mb with QMC-III and WS wave functions, respectively. This is only about 
10% smaller than the lower limit of the corresponding experimental value of 97+10 mb [12]. 
However, since the upper limit of uncertainty in the experimental data is of the order of 20% nl2n . 
this slight under prediction of the total electromagnetic cross sections by our model may not be 
serious. 

In Fig. 4, we show the calculated LMD for the 22 fragment emitted in the CD of 23 O on a Pb 
target at 422 MeV/C beam energy, calculated with QMC-III and WS ground state wave functions. 
The full width at half maximum (FWHM) of the LMD is about 85 MeV/c in both cases. This 



value is more than a factor 2 larger than the FWHM of the LMD of the core fragment observed in 
reactions induced by established one-neutron halo nuclei like n Be and 19 C. Although there are no 
data available for the LMD of the 22 O fragment in the 23 O induced reaction on a heavier target like 
Pb, the FWHM of the LMD of 22 O have been measured in the 23 O induced breakup reaction on 
light targets like carbon [6, Is, 9], where they have been found to lie around 100 MeV/c. Therefore, 
these results almost rule out the development of a one-neutron halo structure in 23 O, even though 
the valence neutron occupies an s-orbit in this nucleus. 



In conclusion, we have studied the structure and the Coulomb dissociation of the neutron drip- 
line nucleus 23 employing the structure wave functions calculated within the latest quark-meson 
coupling model. This is for the first time that a model that invokes the explicitly the quark degree 
of freedom has been used in study of a drip line nucleus. Our calculations confirm a 22 O(0 + )®2si/2 
configuration for the ground state of the 23 nucleus with a one neutron separation energy of 2.86 
MeV in close agreement with the corresponding experimental value. A spin-parity assignment of 
U5/2 to this nucleus is not supported by our model. We predict the development of the neutron 
skin in this nucleus with a thickness of 0.48 fm. 



The valence neutron-core wave function calculated within the QMC model has been used to 
study the Coulomb dissociation of 23 on a Pb target at the beam energy of 422 MeV/c employing 
a theory that requires only the ground state projectile wave function as input and is free from any 
other adjustable parameter. Although at this beam energy the relativistic effects play a role [38|], yet 
a fully quantal relativistic theory of breakup reactions is not yet available. Our theory is essentially 
non-relativistic in nature. Nevertheless, we observe that the existing data on the excitation energy 
spectra and the total electromagnetic one-neutron removal cross section are well reproduced by 
our model. The calculated longitudinal momentum distributions of the 22 fragment are broad 
which effectively excludes the development of a one-neutron halo structure in this nucleus even 
though the valence neutron occupies an s-orbit. 
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